Quantum Transport Theory for Photonic Networks 



Chan U Lei^ and Wei-Min Zhang^'B 

'Department of Physics and Center for Quantum Information Science, 
National Cheng Rung University, Tainan 70101, Taiwan 
(Dated: March 1, 2011, revised) 

In this paper, we develop a quantum transport theory to describe photonic transport in photonic 
networks. The photonic networks concerned in the paper consist of all-optical circuits incorporat- 
ing photonic bandgap waveguides and driven resonators. The photonic transport flowing through 
waveguides are entirely determined from the exact master equation of the driven resonators. The 
master equation of the driven resonators is obtained by explicitly eliminating all the waveguide 
degrees of freedom while the back-reactions between resonators and waveguides are fully taken into 
account. The relations between the driven photonic dynamics and photocurrents are obtained. The 
non-Markovian memory structure and quantum coherence and decoherence effects in photonic trans- 
port are also fully included. This quantum transport theory unifies two fundamental nonequilibrium 
approaches, the Keldysh's nonequilibrium Green function technique and the Feynman- Vernon in- 
fluence functional approach, together to make the investigation of the transient quantum photonic 
transport become more powerful. As an illustration, the theory is applied to the transport phenom- 
ena of a driven nanocavity coupled to two waveguides in photonic crystals. The controllability of 
photonic transport through the driven resonator is demonstrated. 

PACS numbers: 03.65.Yz; 05.60.Gg; 03.65.Db; 42.82.-m; 42.50.Ex 



I. INTRODUCTION 



A photonic network is a communication network in 
which the information is transmitted entirely in terms 
of optical signals With the rapid development of 
nanotcchnology, the network in terms of all-optical cir- 
cuits imbedded in photonic crystals could be a promising 
device to provide a robust interconnect network for op- 
tical communication 0, Q . Photonic crystals are artifi- 
cial materials with periodic refractive index, its photonic 
band gap (PBG) structure together with its characteris- 
tic dispersion properties make the light manipulation and 
transmission much more efficient. For examples, strong 
light confinement can be realized by introducing point 
defect in photonic crystals Q, slow light can be gener- 
ated by waveguide structure with controllable dispersion 
properties, which can be implemented by introducing line 
defect or series of coupled point defects in photonic crys- 
tals d, 0] . Various novel devices have been constructed 
or proposed by the combination of different defects with 
PBG structure, such as optical switches filters [3, 
memory devices |S] and on-chip single photon gun [10|, 
etc. Furthermore, with the new development of the semi- 
conductor nanofabrication technique, many tunable func- 
tional photonic crystal devices have also recently been 
proposed and modeled. Different techniques are utilized 
to characterize various physical quantities in these pho- 
tonic crystal devices. In particular, properties such as the 
resonance frequency of a resonator or the band structure 
of a waveguide can be tuned by changing the refractive 
index of the photonic crystal through thermooptic effect 



electrooptical effect JJl], fluid insertion 13|, or even 
by mechanically changing the structure of photonic crys- 



tals [IJ]. Couplings between different elements in a pho- 
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tonic circuit can be controlled with the help of micro- 
/nanoelectromechanical systems [l^. These dynamical 
tunable devices greatly expand the applications of the 
photonic crystals for photonic integrated circuitry and 
further stimulate the potential application in quantum 
information processing. 

To achieve the goal of quantum information processing 
with all-optical processing, information carriers should 
be individual photons and ultra-fast operations are nec- 
essary. Photonic transmission processes in photonic net- 
works should be able to vary in an extremely quick and 
controllable way. In such a situation, all the photonic 
devices are far away from equilibrium, where the non- 
Markovian memory and quantum coherence and decoher- 
ence dynamics dominate the photonic transport. Thus, 
a fundamental quantum transport theory that can incor- 
porate with the non-Markovian memory and quantum 
coherence dynamics of photons for photonic transmission 
in photonic networks is highly demanded. 

In fact, non-Markovian dynamics in quantum optics 
has been extensively studied for a few-level atom placed 
inside photonic crystals [l^, The typical features 

of the non-Markovian dynamics include the atomic pop- 
ulation trapping (inhibition of spontaneous emission), 
the strong localization of light, the formulation of atom- 
photon bound states and the collective switching behav- 
ior in the vicinity of the PBG [18. - 22j . These features 
can been found by solving exactly the Schrodinger equa- 
tion for the atomic state contained only one photon or 
using a perturbative expansion to the Heisenberg equa- 
tion of motion in powers of the atom-field reservoir cou- 
pling strength. When the number of photons increases or 
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the perturbation dos not work, the problem becomes in- 
tractable and the general non-Markovian dynamics with 
arbitrary number of photons at arbitrary temperature 
for the structured reservoir has not been fully solved. 
Recently, we have utilized the exact master equation of 
a micro/nano cavity coupled to a general thermal reser- 
voir and a structured reservoir in photonic crystals to 
study non-perturbatively various non-Markovian process 
involving arbitrary number of photons at arbitrary tem- 
perature [l^, [13]. However, a quantum photonic trans- 
port theory in photonic networks consisting of all-optical 
circuits incorporating photonic bandgap waveguides and 
driven resonators has not been established in the litera- 
ture. 

A fundamental quantum transport theory should be 
built on a fully nonequilibrium treatment. The mod- 
ern nonequilibrium physics is developed based on the 
Schwinger-Keldysh nonequilibrium Green function tech- 
nique (25I , [26j and the Fey nman- Vernon influence func- 
tional approach [23| . The nonequilibrium Green function 
technique allows a systematic perturbative ^28*, and 
also a non-perturbative 30] study for various nonequi- 
librium phenomena in many-body electronic systems, in 
particular in the steady limit. It has become a very 
powerful tool in the study of steady quantum electron 
transport in mesoscopic physics [3l[. However, such an 
approach has not been utilized to investigate photonic 
transport in all-optical processings in photonic networks. 
Besides, the problem of non-Markopvian memory struc- 
ture and quantum decoherence dynamics has not been 
well explored in terms of the Schwinger-Keldysh nonequi- 
librium Green function technique in the transient quan- 
tum transport. 

On the other hand, the Feynman- Vernon influence 
functional approach [27| has been widely used to study 
dissipation dynamics in quantum tunneling problems 32 1 
and decoherence dynamics in quantum measurement the- 
ory [33} . It is in particular very useful to derive the 
exact master equation for the quantum Brownian mo- 
tion (QBM), achieved by integrating out completely the 
environmental degrees of freedom through the path inte- 
gral where the non-Markovian memory structure and the 
decoherence processes are manifested explicitly [13, [s^j . 
The QBM is modeled as a central harmonic oscillator 
linearly coupled to a set of harmonic oscillators simu- 
lating a thermal bath. Applications of the QBM exact 
master equation cover various topics, such as quantum 
decoherence, quantum-to-classical transition, and quan- 
tum measurement theory, etc. [36l - [38j . However, the 
influence functional approach has not been fully used to 
formulate in a closed form the problem of the quantum 
transport phenomena, although some systematically ex- 
pansion approaches, such as the real-time diagrammatic 
expansion approac h 13911 and the hierarchical equations- 
of- motion apprach (40|, have been developed. Until very 
recently, we have established the exact nonequilibrium 
theory for transient electron transport from the exact 
master equation of nanodevices obtained based on the 



Feynman- Vernon influence functional '41']. 

In this paper, we shall attempt to develop a quantum 
transport theory to depict photonic transport in pho- 
tonic networks, based on the recent development of the 
nonequilibrium quantum theory for the transient elec- 
tron transport dynamics in nanodevices [4l| . The pho- 
tonic network consists of all-optical circuits incorporat- 
ing photonic bandgap waveguides and driven resonators. 
We will focus on the dynamics of photonic transport 
flowing from the resonators into waveguides through the 
controllable photonic dynamics of the driven resonators. 
The photonic dynamics of the driven resonators is deter- 
mined by the exact master equation where the waveg- 
uides are treated as reservoirs. The back-reactions be- 
tween the waveguides and the resonators and thereby the 
non-Markovian dissipation and fluctuation induced from 
the back- reactions are fully taken into account. We fol- 
low the similar procedure of the nonperturbative deriva- 
tion of the exact master equation for nanoelectronics we 
developed recently [13, [i3| to obtain the exact master 
equation for photonic networks. The photocurrents flow- 
ing from the resonators into the waveguides that describe 
the transient photonic transport in the network are then 
obtained directly from the master equation. The full re- 
sponse of the photonic transport under the control of 
external driving fields is explicitly presented. The quan- 
tum kinetic theory based on the Keldysh's nonequilib- 
rium Green function technique is reproduced and also 
generalized within our framework. 

The remainder of the paper is organized as follow. In 
Sec. II, we model photonic networks in terms of a general 
open optical system, and specify it by a photonic crys- 
tal structures in particular. The fundamental Hamilto- 
nian for the photonic network is derived from the first- 
principle. The derivation of the exact master equation 
for the driven resonators based on Feynman- Vernon in- 
fluence functional approach and the calculations of the 
relevant physical observables from the master equation 
are presented in Sec. HI, where the response of the res- 
onator dynamics on the external driving field is also ex- 
plicitly given. In Sec. IV, the photonic transport the- 
ory is established. The photocurrents flowing from the 
resonators into individual waveguides that describe the 
photonic transport in the network are derived from the 
exact master equation, from which the non-Markovian 
memory structure and the coherence and decoherence 
dynamics in the photonic transport can be examined ex- 
plicitly. A comparison of the present theory with the 
electron transport in mesoscopic systems based on the 
Keldysh's nonequilibrium Green function technique is 
also given, from which the generalized lesser (correla- 
tion) Green function which determines completely the 
quantum kinetic theory of photonic networks is explic- 
itly solved. In Sec. IV, we utilize the present theory to 
investigate quantum transport phenomena of a driving 
nanocavity coupled to two waveguides in photonic crys- 
tals, as an illustration. The controllability of photonic 
transport through the driven resonator is analytically 
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FIG. 1: Schematic plots of all-optical circuits incorporat- 
ing photonic bandgap waveguides and resonators, (a) Many 
nanocavities couple to one waveguide at different sites; (b) 
One nanocavity (at the center) couples to many individual 
waveguides. 



and numerically demonstrated. Finally, summary and 
prospective are given in Sec. V. Some detailed deriva- 
tions for the master equation and analytical solutions in 
the weak coupling limit are also given in Appendices. 



II. MODELING THE PHOTONIC NETWORKS 

The photonic network considered in this paper consists 
of photonic band gap structure incorporating with defect 
cavities and waveguides whose spectrums lie within the 
band gap of the photonic crystal, as shown in Fig. [T] 
For the sake of simplicity, we restrict ourselves to pho- 
tonic crystal made up of linear, isotropic and transparent 
medium. Therefore, the photonic network can be char- 
acterized by a real scalar dielectric constant e(r) which 
is explicitly spatial dependent ^2]. To quantize the fields 
in the photonic network in terms of the cavities and the 
waveguides modes, we adopt the system-reservoir quanti- 
zation formalism from [44,] , a generalized canonical quan- 
tization of the electromagnetic field in medium [45| . The 
electric and magnetic fields in the photonic network can 
be expressed in terms of the vector potential A(r, t) and 
a scalar potential (f>{r): 



E 



-V0- 



1 dA 

c~dt 



B = V X A 



(1) 



We work in the Coulomb gauge and in the absence of 
source, which corresponds to the choice of = with the 



transversality condition V • [e(r)A(r,t)] = 0. Then the 
Hamiltonian of the electromagnetic field for the photonic 
network is given by 



-f (V X A{v,t)f 



(2) 



where n(r,i) = e(r)A(r,t)/c'^ is the canonical conjugate 
of the vector potential. 

We shall expand the vector potential in terms of a com- 
plete set of mode functions which are defined by the wave 
equation 



Vx [Vxf;,(r)] 



ffc(r) - , 



(3) 



where the mode functions ffe(r) satisfy the orthonormal- 
ity condition 



j dPve{v)il{v)-h'{v) 



Sk,k' 



(4) 



Utilizing the Feshbach projection technique '46| , we sep- 
arate the whole space of the photonic network by the re- 
gion 5i occupied by the cavity i, the region 5q, occupied 
by the waveguides a, and the region left denoted by 
as the background photonic crystal structure. Because 
of the presence of the PBG structure, the micro/nano 
cavities in the photonic crystal have very high Q factor 
and the loss of the waveguides in the photonic crystal is 
also extremely low, the field in Qf, is then negligible as an 
effect of lossless materials for photonic crystals. Conse- 
quently, the exact eigenmodes of the vector potential in 
the photonic network can be decomposed into the eigen- 
modes in the regions 3^ and 5q,, respectively 



\ " \ " ak 



where Ui(r) is the eigenmode of cavity i, VQfe(r) is the 
the continuum eigenmode of the waveguide a. Both of 
them form a complete and orthonormal basis set in re- 
gions 9i and 9q separately. While a*^ and are the 
expansion coefficients. The eigenmodes of the cavities 
and the waveguides satisfy the following equations with 
the boundary conditions [4J] 



V X [V X u,(r)] - ^^u,(r) 



n X [V X Ui(r)] 



(6a) 
(6b) 



V X [V X v„fc(r)] - 



n X VQfe(r) 



■Vafc(r) = , (6c) 
(6d) 



With the above decomposition of the exact eigen- 
modes, we can expand the vector potential of the electro- 
magnetic field in the photonic network and its canonical 
conjugate in terms of the cavities eigenmodes and the 
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continuous modes of the waveguides 

A(r, t) = Yl i) + II A„ (r, t) 

i a 

n(r, t)^Yl n» (r , + ^ n„ (r, t) 



(7) 



where Ai{r,t),Ili{r,t) and AQ,(r, t), nQ(r, i) are vector 
potential and canonical momentum of the fields in cav- 
ity i and waveguide a, which can be quantized following 
the standard quantum field theory in medium. Without 
the loss of generality and also for simplifying the nota- 
tion, the photonic network being considered consists of N 
single mode cavities and M waveguides, each waveguide 
has a single continuous spectrum (it is indeed straight- 
forward to generalize the result to the case of multi-mode 
cavities). Then we have 



A,(r,t) =c\/-^ [a,u,(r)e-^"'* + atu*(r)e''^"*], (8a) 
Zuii 




n,(r,o 

A„(r,i) 



(8b) 



ie{r) 



2ujk 



[cQfcVafc(r)e" 



-4v:,(r)e'-*], 



cL<fe(i-)e 



(8c) 



(8d) 



where the operators ai,aj and Cak,cl^i; are the creation 
and annihilation operators of the fields in cavity i and 
waveguide a. They obey the following commutation re- 
lations 

[ai,a]] = Sij , [Cafc, cj,,j,,] = Saa'Skk', 

[ai,aj] = [cak^Ca'k'] = [ai.Cak] = [al,CQfc] = 0. (9) 

Substitute the representations (jS]) into the Hamilto- 
nian ^ , one arrives at the following Hamiltonian for the 
photonic network 



Cak 



a k 



^ [ViakO-lcak + ViakaiCak + H.C.] , (10) 



i a.k 



where Viak and Viak are resonant and nonresonant cou- 
pling constant between cavities and waveguides which are 
defined by 



f 

^ [u*(r)xn].[Vxv„fe(r)] , 



^2 



iak 



[Ui{r) X n] • [V X Vafe(r)] 



Here we have ignored the couplings between cavities and 
also between different waveguides because their spatial 
overlaps should be very small. Also, due to the PBG 
structure in photonic crystals, the features of high Q- 
factor micrc/nano cavities and the lossless waveguides 
allows us to neglect the nonresonant terms. Thus, the 
Hamiltonian of the photonic network incorporating M 
photonic bandgap waveguides and N resonators can be 
rewritten simple as 



N 



M 



H hLOia\ai + ^ ^ fl^akclk'^ak 

i—1 a—1 k 

+ H {y^akalco,k + H.c.) , (11) 

i OL^k 

The explicit form of the cavity and waveguide disper- 
sions, LOi and uJak, and the coupling between the cavities 
and waveguides, Viak, depends on the detailed implemen- 
tation of the photonic network in photonic crystals, and 
can be obtained by solving the eigenmode equation ([B]). 
It is worth pointing out that the above derivation of the 
model Hamiltonian can be extended straightforwardly to 
photonic networks imbedded in other metamaterials. 

Furthermore, to control the photonic transport in pho- 
tonic networks, we must incorporate the external driving 
fields applying to the cavities and the waveguides, which 
can be realized by imbedding light emitters such as quan- 
tum dot or nanowire to the cavities and the waveguides 
p. Then the Hamiltonian is modified as 

Hit) =Y,f^^Aa^ + E + 

i i 

Oik ak 

+ hY {V^o.kit)alco.k + V*,kit)cik^'d- (12a) 

iak 

This implementation is very similar to the electron trans- 
port in mesoscopic systems in nanostructures [3l|, [4l| . 
For the electron transport in mesoscopic systems, one 
uses bias and gate voltages applying to the electronic 
leads and the central region to adjust the Fermi surfaces 
of the leads and the energy levels of the central region. 
Here one has to use the external driving fields directly 
acting to the cavities and the waveguides, which leads 
the Hamiltonian to contain explicitly additional terms 
proportional to the external driving fields. 

Mathematically, we can transfer the external driving 
field fak{t) acting on the waveguide a into an equiva- 
lent external field acting on the cavities: fi{t) — >■ fi{t) — 
fi(t) - Viak{t)fak(t)/hujak through a shift of the waveg- 
uide operator Cak Cak = Cq/c + fakit)/hujak- Thus 
the general Hamiltonian for photonic networks can be 
divided into three parts 



(13) 
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where, 

i i 

Heo = ^ fv^akclk^^ak , (14b) 
k 

HTc.it) ^hJ2 {V^ckiMcc,, + V*^,{t)cl,a,) . (14c) 

ik 

The first part Hs (t) is the Hamiltonian of the cavities in 
the photonic network plus the contributions of the driv- 
ing field, fi{t), acting to the cavity i (which also includes 
the driving field fak (t) acting on the waveguides in terms 
of the form of —Viak{t)fak{t)/i^ak, as shown above). The 
second part He is the Hamiltonian of the waveguides 
in the photonic network. The third part Hxit) is the 
coupling between the cavities and the waveguides, which 
can be time dependent in general through the controls 
of the coupling between different elements, by means of 
dynamically tuning the photonic structure ^, J5J . With 
the help of the general Hamiltonian p^ . we can model 
various photonic networks. For the sake of simplicity, we 
take h=l hereafter. 



III. EXACT MASTER EQUATION FOR 
DRIVEN RESONATORS 



Experimentally, the photonic dynamics of resonators 
and the photonic transport between waveguides and res- 
onators can be controlled by the external driving field. 
When a resonator has exchanges of particles, energy and 
information with the surroundings, it becomes a typical 
open system. For an open quantum system, its dynam- 
ics cannot be properly described by Schrodinger's wave 
equation. It should be described by the master equa- 
tion of the reduced density matrix [i^. The reduced 
density matrix of the resonators, denoted by p{t), fully 
depicts the dynamics of photonic coherence of the driven 
resonators coupled with many waveguides. Any other 
physical observable is simply given by {0[t)) — {Op{t)). 

The reduced density matrix describing the photonic 
quantum state of the driven resonators is defined by trac- 
ing over all of the waveguide degrees of freedom from the 
total system (the driven resonators plus waveguides): 



p{t) = trs [ptot (i)] 



(15) 



where ptot(i) is the density matrix of the total system. 
The total density matrix follows formally the evolution 
equation: 

Ptot(<) = U{tM)ptot{to)U\t,tQ) , (16) 
where U{t,to) is the evolution operator, 

U{t,to) = Te^p{-i H{T)dT}, (17) 



in which T is the time-ordering operator, and H{t) is 
the Hamiltonian of the total system given by Eq. (|13|) . 
In the following, we will derive the exact master equation 
for the reduced density matrix by explicitly integrating 
out all of the waveguide degrees of freedom based on the 
Feynman- Vernon influence functional 27], similar to the 
derivation of the master equation for electron coherence 

tiamics and electron transport in nanostructures (4l| - 
and also for a coupled cavity fields in a general bath 
. The difference is that here the system contains ex- 
ternal driving fields and the reservoir is at arbitrary finite 
temperature initially, which makes the derivation much 
more complicated. 

As usual, we assume that there is no initial correla- 
tion between the resonators and waveguides [s^ , namely 
Ptot(to) = p{tQ)® PE(to), and the waveguides are initially 



PEito) 



in the equilibrium state, i.e. 
with Pa = ^/iksTa) and Tq, is the initial temperature 
of the waveguide a. Then the reduced density matrix at 
arbitrary later time t can be expressed in the coherent 
state representation as [i^ \^ : 

{af\p{t)\a'j) 

= J dfi{ao)d^{ctQ){ao\p{to)\aQ)J{a*j:, a'f,t\ao, ct'l, to) 

(18) 

with the vector a = (ai,a2,---) and \ck) — |q;i)|q!2)--- is 
a unnormalized multi-mode coherent states, i.e. ai\ci) = 
ai\a.) and {ot\a') = exp(^-a*a-) = exp(Q;'''a'), while 

dfj,{oi) ~ ^"^'^" e~l" is the integral measure of the 
Bergmann complex space [i^]. The propagating func- 



tion of the reduced density matrix in Eq. ()18p is given in 
terms of path integrals: 

J{ct*j,a'j,t\ao,a'l,to) 

(19) 

with the integral boundary conditions: Q;(io) = cto, 
a*{t) — cx*^, a'*{to) = (x'q, and OL'{t) = a'f, where 
S'5[a*,a] is the action of the driven resonators in the 
coherent state representation: 



Ss[a\a\^-^-[oi\a.{t) + a.\ta)oLQ]+ j dr|^ 



, da. 
d7 



dOL^ 



(Q;^a;a + at/ + /ta)|, (20) 



in which the frequency matrix u; = {w^} (here the off- 
diagonal matrix elements representing the possible cou- 
plings between different cavities are included for the gen- 
erality of the formulation) and the driving field vector 
/(r) = {fiM}- While ^[a*a;a'*a'] is the influence 
functional [27] obtained after integrated out all of the 
waveguide degrees of freedom l47j : 
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^[a*a;a'*a'] =exp{^(- f' dr r dT'a\T)ga.{r,T')niT') - t dr r dr' a.'\T')^c.{T' ,t)cx' [t) 

Q, t-O J to to J tQ 

+ f dT f dT'cy'\T)g^{T,T')cx{T')- f dr rdr'[at(r)-a't(r)]g„(r,r')[a(T')-a'(r')])}. (21) 



The time-correlation functions in Eq. (j21p are given by: 
g»,(r,r') = ^T^,afc(r)F/„,(T')e— (^-^') , (22a) 

A; 

gay(r,r') = ^y,„fc(T)V/„,(r')n„(u;„fe)e-^--(^-^'), 

(22b) 



which depict the time correlations of photons in the 
waveguides through the resonators, and ^^(wqj.) = 
l/(e^""°'' — 1) is the initial thermal photonic distribu- 
tion function in the waveguide a at the time t^. The 
influence functional, Eq. ((2T|) . contains all the back reac- 
tions from waveguides to the driven resonators through 
the photonic transfer between the resonators and waveg- 
uides. If the coupling constants do not explicitly de- 
pend on the time, then introducing the spectral density 
Jq(w) = '2.TTY,k^iakV*^f,5{uj - ujak), wc Can simply ex- 
press the time-correlation function in terms of the spec- 
tral density as follows 



gmi(T,r') =y ^Ja(cj)e" 



o{t-t') 



laijir, t') 



duj 
2^ 



-iuj[r — r') 



(23a) 
(23b) 



As shown in Eqs. (PT|) and (IT^ . after integrated out all 
of the waveguide degrees of freedom, the effective action 
in the propagating function remains in a quadratic form. 
Therefore, the path integral in the propagating function 
can be carried out exactly through the stationary path 
method. The resulting propagating function, Eq. (|19p. 
becomes (see a detailed derivation given in Appendix A) 

J{a*f., a'f, t\an, a'o, to) = M^) cxp{(Q:^ - y\t))J i{t)oLa 
+ a'lj\{t){a'j - y{t)) + a\j2{t)a'f + a'J J3(t)ao 
- a^yit) - y\t)cx'f - yHt)w{t)y{t)}, (24) 

in which A{t) = det[w{t)], Ji{t) = w{t)u{t,to), Jiit) = 
l-w{t), Jsit) = l-u^t,to)w{t)u{t,to) andu;(t) = [1 + 
v{t,t)]~^ = w^t). The function m(t, to), u{T,t), v{T,t) 
are TV by iV matrices, and y{t) is a 1 x matrix, where 
N is the number of resonator modes. These functions 



I 

obey the equations of motion 

du{T, to) 
dr 

du{T, t) 



dT 

dv{T, t) 
dT 

dT 



iu:u{T, t)- I g(T, t')w(t', t)dT' = 



(25a) 



(25b) 



iu3v{T,t)+ I g{T,T')v{T',t)dT' 
'to 

rt 

1{t,t')u{t' ,t)dT' , (25c) 

/to 

li^y{T)+ f g{T,T')y{T')dT' =-lf{T) (25d) 
Jtn 



subjected to the initial conditions u{t,t) = 1, tt(io,to) = 
1, v{tQ, t) — and y{to) ~ with tg < t,t' < t. Here we 
have also defined g(T, r') = Xla ''"') ^^'^ si''','''') = 
J2aSa{T,T'). It is not too difficult to find that 



(26a) 
(26b) 



u{T,t) = u\t,T) , 

yir) = -t r u{T,T')f{T')dT' , 
Jto 

v{T,t)= f dT' f dT"u{T,T')g{T',T")u{T",t) . (26c 
Jto Jto 



As we see, once Eq. (|25ap is solved for it(T, tg), the dy- 
namics of the driven resonators can be completely deter- 
mined. 

Substituting Eq. ([M]) into Eq. ([TH]), taking a time 
derivative to the reduced density matrix, and then using 
the relations [i^: ai\a.) = a^ja) and -^\a.) = a|^|Q;), we 
obtain a time convolutionless but exact master equation 
for the driving resonator system coupled to waveguides: 



dpjt) 
dt 



+ 7»j(t)[2ajp(t)aJ - p(t)a|aj - ajajp(t)] 

(27) 

where i?eff(t) = (*)«!%■ + EJ/((0«1 + /'*(*)«.] 

is the effective Hamiltonian of the driven resonators. 
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The renormahzed frequencies uj'ijit) and the renormal- 
ized driving field //(i) result from the back- reaction ef- 
fect of the waveguides coupled to the resonators. The 
initial temperature effects are fully incorporated into the 
temperature-dependent noise coefficient 'jit). The time 
dependent dissipation coefficients 7(t) together with the 
time-dependent noise coefficients 'jit) describe the non- 
Markovian dissipation and decoherence dynamics of the 
resonators due to the interaction between the resonators 
and waveguides. Here these coefficients also describe the 
photon transport from the resonators into the waveguides 
under the control of external driving fields, as we will 
show in the next section. All these time-dependent coef- 
ficients in Eq. ^27\i are given explicitly as follow: 



a 



(28a) 

(28b) 
(28c) 
(28d) 



with 



Ka{t) = / dTga{t,T)u{T,to)u ^{t,to), (29a) 

Kit) = / dT[gait,T)v{T,t) -ga{t,T)u{T,t)] 
Jtn 



Ka{t)v{t,t) , 



(29b) 



f^{t)^iK4t)y{t)-i f dTga.it, T)y{T) (29c) 

as functions of u{t,to), v{t,t) and y{t) that are deter- 
mined by the intcgrodifferential equations of Eq. 
We should point out that although it has a time convolu- 
tionless form, the master equation is exact and the non- 
Markovian memory effect between the resonators and 
waveguides is fully embodied non-perturbatively in the 
integral kernels involving the non-local time-correlation 
functions of the waveguides, g(T, r') and g(T, r'), in 
Eq. (P5|l . We should also point out that /„(*) is a shift 
to the external driving field f{t), which comes from the 
back-reaction of the waveguide a to the resonators. In 
other words, /^(i) is a feed-back effect to the exter- 
nal driving field. We may call the function y{t) (see 
Eq. (j26bp ') as the driving-induced field. 

Notice that using the Feynman- Vernon influence func- 
tional approach to derive the exact master equation was 
carried out early for quantum Brown motion (QBM) 
[35I [50| . The exact master equation can also be obtained 
by the trace-over-bath on total-space Wigner-function 
method l5l[453ll or the stochastic diffusion Schrodinger 
equation j54{ . All these previous works deal mainly with 
a single harmonic oscillator coupled to a thermal bath 



without external driving fields. The extension of the 
exact master equation to the systems of two entangled 
optical modes or two entangled harmonic oscillators has 
only been worked out very recently [13, . The ex- 

act master equation of Eq. ((27)) obtained here is indeed 
the most general one for the system containing arbitrary 
number of entangled modes coupled to arbitrary number 
of photonic reservoirs with arbitrary spectral density at 
arbitrary initial temperatures under arbitrary number of 
external driving fields. It allows to investigate various 
general and exact non-Markovian dynamics in photonic 
systems. 

On the other hands, one may obtain a time- 
convolutionless master equation from perturbation ex- 
pansion up to the second order in terms of the cou- 
plin g c onstant between the system and the reservoir 
[36l . |48|. Such a time-convolutionless master equation 
is valid only in the weak coupling regime where the non- 
Markovian effect is also very weak, except for some struc- 
tured reservoirs Also, the time-convolutionless per- 
turbation master equation up to the second order in the 
system-reservoir coupling is indeed the same as the Born- 
Markov master equation without taking the long time 
limit [23, . An exception of the time convolutionless 
exact master equation not for pure optical systems is the 
master equation of a two-level atom with one single pho- 
ton at zero temperature (stI . [ssj . However, when the 
number of photons increases, the problem becomes in- 
tractable and the general non-Markovian dynamics with 
arbitrary number of photons at arbitrary temperature 
has not been well explored in quantum optics. Different 
from the study on the non-Markovian atomic dynamics 
involving only one single photon in quantum optics, the 
exact master equation of Eq. ([27| is capable to investi- 
gate the general and the exact non-Markovian process 
involving arbitrary number of photons of the reservoir 
with arbitrary spectral density at arbitrary temperature 
in all-optical circuits. Applications to such exact non- 
Markovian dynamics of a micro/nano cavity coupled to 
a general thermal reservoir and a structured reservoir in 
photonic crystals have just been worked out very recently 

MM. 



In fact, the master equation ([27)) determines all the 
nonequilibrium dynamics of the driven resonators and 
the photonic transport between the resonators and the 
waveguides in photonic networks. The aim of the present 
paper is to use the exact master equation derived above 
to establish a quantum transport theory that generalizes 
the quantum transport theory based on the Keldysh's 
nonequilibrium Green function technique. We will show 
in the next section how the photocurrent passing through 
the resonators into waveguides can be obtained directly 
from the master equation, and how the quantum trans- 
port theory based on the Keldysh's nonequilibrium Green 
function technique can be easily reproduced and general- 
ized. Before end this section, we shall calculate some im- 
portant physical observables describing the photonic dy- 
namics in the driven resonators. One of them is the time 
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evolution of the photonic resonator field, and another one 
is the single particle reduced density matrix which char- 
acterizes photonic intensity of the resonators. These are 
two main quantities related to the non-Markovian trans- 
port dynamics in photonic systems. 

The photonic field of the resonators is determined by 
{ai{t)) ~ tis[aip{t)]. With the help of the master equa- 
tion, it obeys the equation of motion: 



dt 



{a^ =Y,[u{tM)u-\tM)h{aj{t))-ifl{t) . (30) 



Its solution is: 

{a,{t)) =J2^^Jit,to){aJ{to))+y^{t) . (31) 



In other words, the photonic field of the driven resonators 
is a combination of the photonic propagating of the initial 
field u{t,tQ){aj{to)) and the driving-induced field y{t) 
(see Eq. (j26b|) '). This solution describes the driving- field- 
induced photonic coherence in resonators. Similarly, the 
equation of motion for the single particle reduced density 



'-^''^^-'^ — ti g[a'^aip{t)], can also 
be found from the exact master equation: 



matrix of the resonators, p\j {t) --s l-j 



m 

-in{t){a]{t))+i{a.,{t))r^\t) . (32) 

The solution of the above equation of motion can also 
be obtained explicitly in terms of the functions u{t,tQ), 
v{t,t) and the driving-induced field y{t): 

mn 

+ {a,(t))v*{t)+v,{t){a]{t)) -y,it)y*{t). (33) 

The first term in the above solution comes from the initial 
photonic distribution in the resonators. The second term 
is the effect induced by the thermal fluctuation in the 
waveguides (see the solution of v{t, t) by Eq. (I26cp '). The 
remaining terms are the contribution coming from the 
driving fields. Besides, the ith diagonal element of the 
single particle reduced density matrix gives the average 
photon number of the mode i in the resonators, i.e. 



p'i\t)^tT,[a^^a,pit)]^n,{t), 



(34) 



which measures the photonic intensity of the ith mode in 
the resonators. The photonic coherence and the photonic 
intensity of the driven resonators, given respectively by 
Eqs. (|3T1) and (p4|). play the important role for the pho- 
tonic transport in photonic networks, as we shall discuss 
in the next section. 



IV. PHOTONIC TRANSPORT IN PHOTONIC 
NETWORKS 

In this section, we are going to derive the exact pho- 
tocurrent flowing from the resonators into each waveg- 
uide. Connection between the photocurrent and the mas- 
ter equation of the reduced density matrix is given, which 
explicitly demonstrates the intimacy between quantum 
coherence and quantum transport in nonequilibrium dy- 
namics. Also, the relation between the transport the- 
ory obtained here with the transport theory based on 
the Keldysh's nonequilibrium Green function technique 
is discussed and the powerfulness of the present theory is 
given. 



A. Photocurrent in each waveguide channel 

To flnd the photocurrent flowing from the resonators to 
each waveguide using the master equation, it will be more 
convenient to reexpress the master equation, Eq. (j27p . as 
follows: 



dpjt) 
dt 



= -i[Hs{t),p{t)]+Y,[^+it)+C-itMt) . (35) 



Here Hs{t) is the original Hamiltonian of Eq. (|14ap for 
the driven resonators, /^^(i) and C'^it) are superoper- 
ators acting on the reduced density matrix, induced by 
the coupling to the waveguides. They are given by: 



- K.aij{t)a\ajp{t) - ijai{t)a\p{t) + H.c.} , (36a) 
C WpW = ^{^ctj{t)[a\p{t)aj - a\ajp{t)] 

+ Ka,,j{t)ajpit)al + if^,{t)p{t)a\ + H.c.} (36b) 

The superoperators C'^{t) and C'^it) are intimately re- 
lated to the photocurrent through the waveguide a, as 
we will see next. 

The photocurrent flowing from the resonators into the 
waveguide a is dcflncd in the Heisenberg picture as: 

^ iY}y^o.J.{t){a\{t)co.u{t)) - V:^^{t){a\{t)co.u{t))] , 

ki 

(37) 

where Na — '^ak'^ak is the photon number operator 
of the waveguide a. Furthermore, consider the single 
particle reduced density matrix in the Heisenberg picture: 
p[]\t) = tYs[aUip{t)] = {aUt)ai{t)), we find that 



dp(^^t) 
Jt 



i[iv,p^'Ht)]+S{t)~J2x^{t) , (38) 
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where Ta{t) and S{t) are the current matrix of the 
waveguide a and the source matrix of the driven res- 
onators: 



S,,{t)^^{a,{t))fnt)-Km{aUt)) , 



(39a) 



(39b) 



Comparing Eqs. (|37| and (j39bp . one can see that the 
trace of the current matrix is just the photocurrent flow- 
ing from the resonators into the waveguide a: Ia(t) = 
TT[Xa{t)]. Note that Tr is the trace over the N x N ma- 
trix of the N single modes in the resonators, while tr^ and 
tr^ used before denote the traces over all the quantum 
states of the resonators and the waveguides, respectively. 

On the other hand, the equation of motion for the sin- 
gle particle reduced density matrix can also be obtained 
directly from the master equation, Eq. (|35l) . The result 
is 



z[a;, p(i) (<)],, +Sy(i) 
Y,^r,[a]aACl{t) + C-{t)]p{t)] . (40) 



Comparing Eqs. (|38p and (|40l) for the single particle re- 
duced density matrix, with the help of Eqs. p6p and 
we obtain the explicit formula for the photocurrent 
matrix: 

X„(i)= [ dT{g„(i,T)p(i)(T,i)-g„(i,r)w(r,t) + H.c.} 

(41) 

and p^^^Tjt) is the generalized correlation function 
which is given explicitly by: 

inn 

+ yi{'T)y*j{t) + ^[uim.{T, ta){a,n{tQ))y*{t) 

ni 

+ y.M(«L(io)X,(i,io)] . (42) 

Then the photocurrent flowing into the waveguide a is 
simply given by: 

J„(t) = 2Re I dTTr[g„(t, t)p'^^\t, t) - ga.{t, t)u(t, t)]. 

(43) 

It shows that the photocurrent is completely determined 
by the time-correlation functions ga(i,T) and gQ(i,T) 
of the waveguides and the propagating and correla- 
tion functions u{T,t) and p^^\T,t) of the driven res- 
onators. The time-correlation functions characterize the 
non-Markovian memory structure of photon transfer be- 
tween resonators and waveguides, while the propagating 
and correlation functions depict completely the photon 
coherence and photon correlation of the resonators un- 
der the control of external driving fields. 



In fact, Eq. (j38|) is a generalized quantum continuous 
equation. By tracing over the equation of motion for the 
single particle reduced density matrix, we have 



(44) 



Here N{t) — Tr[p'^^^(t)] is the total photon number in the 
resonators, S{t) = Tr[<S(i)] is the source coming from the 
driving fields, and Ia{t) = TT[Xa{t)] is the photocurrent 
flowing into the waveguides a. Eq. tells that the 
increase of the photon number in the resonators equals to 
the receiving photons from the driving field subtracting 
the lessening photons flowing into waveguides. 



B. Relations to the Keldysh's nonequilibrium 
Green function technique 

As one seen from Eq. the photocurrent is com- 

pletely determined by the time-correlation functions 
ga{t,T) and ga{t,T) of the waveguides plus the prop- 
agating and correlation functions u{T,t) and p^"'"^ (r, t) 
of the driven resonators. In fact, we have shown |41| 
that the functions tx(r, to), u{T,t), p^^\T,t) are related 
to the retarded, advanced and lesser Green functions of 
the resonators in the Keldysh's nonequilibrium formulism 



(45a) 

u^J{h,t2) ^ e{t2 - ti){[a,{ti),a]{t2)]) = -iG^j{ti,t2), 

(45b) 

p'i]\h,t2) = (a^(^2)a^(^l)) = -^G<(^l,^2). 



(45c) 



While, the time-correlation functions ga{t,T) and 
ga{t, t) correspond to the retarded and lesser self-energy 
functions arose from the couplings between the res- 
onators and the waveguides: 



laij{ti,t2) = iS^jj-(ti,t2), 
la.^J{t^M)=-i^U{tl.t2)■ 



(46a) 
(46b) 



The explicit form of these self-energy functions is given 
by Eq. (1^^ or (|23l) which is rather simple. In 
the Keldysh's nonequilibrium Green function technique, 
quantum transport theory is completely determined by 
the retarded, advanced and lesser Green functions. 

Explicitly, Eq. p5a|) for u{T,to) obtained in the last 
section can be rewritten as: 



+ f Y.^{T,T')G'-{T'M)dT' 

Jtn 



(47) 



which is just the standard Dyson equation for the re- 
tarded Green function. The advanced Green func- 
tion obeys the relation: G°'{ti,t2) = [G'^(i2, ^i)]^, see 
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Eq. (I26ap . The central and also the most difficult part in 
the Keldysh's nonequilibrium Green function technique 
is the calculation of the lesser Green function G^(r, t). 
The lesser Green function G^{ti, t2) fully determines the 
quantum kinetic theory of nonequilibrium system. From 
Eq. (|42|) . we have obtained already the exact analytical 
solution of the lesser Green function: 

G<{T,t)^G''{T,to)G<{hM)G^{to,t) + iy{T)yHt) 

+ r dn I dr2 G''(T,Ti)S<(ri,T2)G"(r2,t), 

(48) 

where G^{to,to) = ip^-^\to) is the initial photon distri- 
bution in the resonators, (a(to)) is the initial resonator 
fields, and y{t) is the driving-induced resonator field 
given by Eq. (j26bp . 

Comparing with the electron transport in mesoscopic 
systems, one always has (a''' (to)) = (a(to)) ~ for elec- 
trons. Also, one usually ignores the first term in the 
standard Green function calculation by taking to ^ —oo 
which loses the information of the initial state depen- 
dence in quantum transport, an important effect on non- 
Markovian memory dynamics. Besides, the external driv- 
ing fields applied to the leads and gates in mesoscopic sys- 
tems are embedded into the spectral densities or the en- 
ergy levels in the central region so that no extra driving- 
induced field is produced, namely y{t) = 0. Thus, the re- 
sulting lesser Green function obtained in the mesoscopic 
electron transport contains only the last term in Eq. 
[Slj . In other words, Eq. gives the exact and gen- 
eral solution for the lesser Green function in photonic 
systems. With the above relations and solutions, the 
photocurrent, Eq. ([iS]). can be re-expressed as: 

Ut) = 2Re / dTTr[S^(t, r)G<(r, t) + S<(t, r)G"(r, t)]. 

Jto 

This reproduces the standard transport current in the 
Keldysh's nonequilibrium Green function technique that 
has been widely used in the investigation of various elec- 
tron transport phenomena in mesoscopic systems, al- 
though most of the previous works used only the special 
solution of the lesser Green function, as we just discussed 
above. 

In conclusion, a full quantum transport theory for pho- 
tonic dynamics in photonic network has been established 
based on the Feynman- Vernon influence functional ap- 
proach. In the literature, the investigation of quan- 
tum transport used mainly the Keldysh's nonequilibrium 
Green function technique. The Keldysh's nonequilibrium 
Green function technique has the advantage of treating 
the lesser Green function in a complicated system by the 
assumption of adiabatically switching on the many-body 
correlations. This allows one to trace back the initial 
time to — — oo, which provides a great simplification for 



practical evaluation but meantime it lacks a proper treat- 
ment of transient phenomena. The Feynman- Vernon in- 
fiuence functional aims to address dissipative dynamics 
of an open system in terms of the reduced density matrix 
(an arbitrary quantum state). The master equation de- 
rived from the influence functional explicitly determines, 
by definition, the temporal evolution of an initially pre- 
pared state. For the photonic networks considered in this 
paper, the many-photon correlations are less important 
and can be ignored so that we are able to obtain the exact 
master equation where all the non-Markovian memory 
effects are encoded into the time-dependent coefficients 
in the master equation. It turns out that these time- 
dependent coefficients are determined indeed by the re- 
tarded and less Green functions in Keldysh's formulism, 
as we have just shown. Therefore, all the advantages of 
the nonequilibrium Green function technique are main- 
tained in our theory but the difficulty in addressing 
the transient dynamics in Green function technique is 
avoided in terms of the master equation. As a result, we 
unify two fundamental nonequilibrium approaches, the 
Schwinger-Keldysh nonequilibrium Green function tech- 
nique and the Feynman- Vernon influence functional ap- 
proach, together to make our theory more powerful in the 
study of the transient transport phenomena in nonequi- 
librium photonic systems. 



V. ANALYTICAL AND NUMERICAL 
ILLUSTRATION 

A. Waveguide as a tight-binding model 

In this section, we apply the theory developed in the 
previous sections to a photonic circuit in photonic crys- 
tals as an illustration. One of the physical realization for 
the photonic network described in Sec. II is the system of 
one or a few nanocavities coupled to many waveguides in 
the photonic crystals, as schematically shown by Fig. [1] 
Here a nanocavity can be considered as apoint defect cre- 
ated in photonic crystals as a resonator [j] . Its frequency 
can easily be tuned to any value within the band gap by 
changing the size or the shape of the defect. While, a 
waveguide in photonic crystals can also be considered as 
a series of coupled point defects in which photon propa- 
gates due to the coupling of the adjacent defects 
By changing the modes of the resonators and the cou- 
pling configuration through various techniques (lTl - [l3 |. 
the transmission properties of the waveguide can also be 
manipulated. In principle, one can solve Eq. ^ to ob- 
tain the dispersion relations of the waveguides and the 
couplings between the cavities and waveguides 0, [l0| . 
Here we may treat the waveguide as a tight-binding 
model, namely the Hamiltonian of the waveguide and 
the coupling Hamiltonian between the nanocavities and 
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the waveguide can be expressed explicitly as |24 |: 

Heu = ^ WQc]j„Can " ^ Ca(cJ,„CQ„+i + H.C.) 



i?Ta = ^mn. (ajccm, + H.C. 



(50a) 
(50b) 



where is the hopping rate between adjacent resonator 
modes within the waveguide a, and it is experimen- 
tally turnable. ^iam is the coupling constant of the zth 
nanocavity coupled to the nth resonator in the waveguide 
a, which is also controllable by changing the geometrical 
parameters of the defect cavity and the distance between 
the cavities and the waveguide !5^. 

Furthermore, making the Fourier transform for semi- 
infinite waveguide, 



(51) 



Cafc = Y - ^ sin(nA:)cQ 



we can transform the Hamiltonian of Eq. (1501) from the 
spatial space into the wavevector space with the result: 

Hec = '^iOakcl^Cak , (52a) 

fe 

Hra = {V^a.kalcak + H.C.) , (52b) 

k 

where < fc < tt, and uJak and Viak are given by: 



UJak =^a- 2^Q COs(fc), Viak 



^iam sin(nifc), (53) 



and cl^^. , Cak are the creation and annihilation operators 
of the corresponding Bloch modes in the waveguide a. 

More specifically. Fig. [ija) consists of N nanocavities 
coupled to one waveguide at different sites {rii}. Then 
the total Hamiltonian of the system can be rewritten ex- 
plicitly as: 



N 



Ha{t) (w.ala, -t- ma] + f*{t)a,) 

=1 

ujkclck + Y {Vtkalck + H.C.) , (54) 



i=l 



ik 



where cok = ujq- 2^o cos(fc), Vik = y f sin(nifc). The 

network presented in Fig. Wijo) shows a nanocavity (at 
the center) coupled to M individual waveguides. The 
corresponding Hamiltonian is 

Hbit)=LUca^a + f{t)a^+r{t)a 

+ ^a;Qfccj^^Cafc +^ (Vafca^Cafc +H.C.) , (55) 



ak 



ak 



where the band structure of the waveguides and the cou- 
pling constants are given by Eq. (j53p with = 1. These 



CO, 



T 



CO., 
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FIG. 2: A nanocavity with frequency lo^. couple to two coupled 
resonator optical waveguides (CROWs). 



physical realizations specify the model Hamiltonian of 
Eq. (HSl). 

As an illustration, we consider here simply a driven 
nanocavity coupled to two coupled resonator optical 
waveguides (CROWs) in photonic crystals, as plotted in 
Fig. [21 and then calculate analytically and numerically 
the photonic transport phenomena. The Hamiltonian of 
this simple photonic circuit is : 

R = uj^a^a + (i;oe-"^''*at + E^e^'^'^U) 

2 2 

(56) 



a=l k 



Q=l k 



where i?o is the strength of the external driving fields in 
frequency a;^. The frequency uj^k and coupling constant 
Yak are given by Eq. (|53l) with rii = \. The corresponding 
spectral densities are given by Ja{i^) = '^''^ga{i^)\Va{i^)\'^ 
where (?(w) is the density of state of the waveguide a and 
Va{u)) is the coupling between the cavity and the waveg- 
uide a. They can be calculated directly from Eq. ([55]) : 



dk _ 1 

1 /Cc 



(57a) 



:(..) = -=(|^)V4e^F^^, (57b) 



with LUa — 2^Q < uj < uja + 2^q . Then the spectral density 
can be explicitly written as 



J (Uj)^} vlV'^^a- {^~^c.y , |W - UJa\ < 2^ 
\ , \uJ~UJa\> 2^a 



(58) 



with r]a = ^ca/^a- In practical, ^ ^a, namely the 
waveguide has a very narrow band. 
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B. Analytical solutions in the weak coupling regime 

First, we shall discuss the cavity photonic dynamics. 
Consider the case where the nanocavity is initially empty: 
(a(io)) = and n{to) = 0. The cavity field and the cavity 
photon nmubcr at a later time, i.e. Eqs. (|3T|) and (f33| . 
become 



Ht)) = y{t) , 

n{t) = {a\t)a{t)) ^v{t,t) + \y{t)\' 



(59a) 
(59b) 



Here y{t) is the driving-induced field given by Eq. (|26bl) 
and v{t,t) is the photon correlation due to the thermal 
fiuctuation in the waveguides, given by Eq. (|26cl) . Only 
at zero temperature, we have v{t, t) = 0. Then the cavity 
photon number equals to the absolute square of the cavity 
field, i.e. n{t) = \{a{t))\'^. 

To see the dynamics of cavity field controlled through 
the driving field, we shall consider first the cavity decou- 
pled to waveguides, i.e. ^ic = ^2c = 0. In this situation, 
the photonic propagating function in the cavity is simply 
given by u{t, to) = e"*""^'*"*"' and the photon correlation 
function v{t,t) = 0. Then, the cavity field and cavity 
photon number equal to y{t) and |y(i)p, respectively, 
with 



(60) 



It shows that the cavity field is a coherent superposition 
of the cavity mode with the external driving field. If 
the driving field frequency differs from the cavity mode, 
the interference between the cavity mode and the driving 
field causing the photons jump in and out of the cavity 
with the Rabi frequency {ujd — a;c)/2. When the driving 
field is in resonance with the cavity mode, the driving 
field would be completely absorbed into the cavity, and 
the cavity field amplitude increases linearly in time with- 
out an interference oscillation. 

When the cavity couples to waveguides, generally it 
is not easy to find the analytical solution for u{t,tQ) and 
v(t, t). One has to solve Eqs. (j25bl) and (I26c[) numerically 
to understand the photonic dynamics of the driven cavity. 
However, in the weak coupling regime, the memory effect 
is negligible so that the Born-Markov approximation can 
be applied [23| . Then the photonic propagating function 
and the correlation function in weak coupling limit can 
be approximated (see Appendix B) by 



j,-(icj^ + fc)(t-io) 



u{t,to) 



_ -2K(t-io)l 



(61a) 
(61b) 



where uj'^ = i^c + ^'^a is a renormalized cavity fre- 
quency with frequency shift SuJa ^ V[f f^^^]- The 
damping rate k = with — Ja{uJc)/2 and the 



average photon number of the two waveguides n(uj[, , T) = 

Ea>^a(wO"a(Wc)/2K. 

On the other hand, to solve the driving-induced field 
y{t) in the BM limit, special care needs to be taken since 
the driving field has its own character frequency which 
is usually different from the cavity mode frequency Wc. 
Therefore, instead of applying directly the BM solution 
of Eq. (|61ap into Eq. (j26bp , we have to use different char- 
acter frequencies ujc and cod for the homogeneous and in- 
homogeneous solutions of Eq. (I25dp to find the BM limit 
of the driving- induced field y{t). The detailed derivation 
is also given in Appendix B. The result is nontrivial: 



where 



= tan 



with the driving field induced frequency shift Suia — 

-P^^tw^-, « = ELi^" with = Uud)l2, lo- 
cated at the driving frequency rather than the cavity fre- 
quency. Correspondingly, the cavity field and the cavity 
photon number of Eq. ([55]) in the BM limit become 



i{t)) 



-icjdit—to) 



-{iuj'^+K){t~to)] 



(63a) 



_ p-2K(t-to)l 



{uJd - Wc)2 + 

_ 2e-(*-*o) COS [{LJd~io'J{t^to)] 



(63b) 



As one see, the cavity field is a coherent superposition of 
the external driving field and the damped cavity mode, 
where the damping comes from the coupling of the cavity 
with waveguides that induces photon dissipation from the 
cavity into waveguides. This damping (dissipation) effect 
leads cavity mode to vanish in the steady limit. Only the 
driving field is remained in the cavity with a modified 
(amplified) field amplitude Eq — ^° and a 

phase shift (j), as a feed-back effect from the coupling of 
the cavity with waveguides. The photon number in the 
cavity is then a combination of the driving-induced field 
plus a background noise from the thermal fluctuation of 
waveguides n{uj'^,T)]. These analytical results help us 
to understand the driven cavity dynamics in the strong 
coupling regime. The corresponding numerical solution 
will be presented later. 

Furthermore, the transport phenomena of photons 
flowing into the waveguides can be specified by the pho- 
tocurrent which can be simplified as well in the weak 
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coupling limit: 



{K,Ql 



g-2K(t-io) 



(1 + -ii)e-«(*-*o) COS [{u:d - u:',){t - to)] 



in[(a;d-a;^)(t-to)]|- 
(64) 



The first term in Eq. is the contribution of the 

waveguide thermal fluctuation, which would be vanish 
when the cavity is equilibrated with the waveguides. 
The remainders come from the response of the waveg- 
uides to the external driving field through the cavity, 
in which the first decay term is due to the dissipa- 
tion of the cavity mode, the second and third decay 
terms are the thermal-fluctuation-induced decoherence 
(noise) of the cavity field. These three contributions to- 
gether with the thermal-fluctuation-induced current (the 
first term) will vanish in the steady limit. The time- 
independent term in Eq. (|64p is the steady photocurrent, 
which is determined by the cavity mode, the driving field 
frequency and the spectrum of the waveguide. When 
the driving field frequency lies outside the band of the 
waveguide a, the steady photoncurrent vanishes because 
K,a = Jai^d)/^ = 0. On the other hand, the photocur- 
rent flowing into the waveguide a becomes maximal when 
the driving field frequency lies at the band center of the 
waveguide a, i.e. uJd = ^a, see Eq. ([55]) . This result 
shows that the nontrivial Born-Markovian solution in the 
presence of external driving field gives already a clear pic- 
ture on the controllability of photonic transport in the 
photonic circuit through the external driving frequency. 



C. Exact numerical solutions in both the weak and 
strong coupling regimes 

Now, we turn to the exact numerical calculation for 
arbitrary coupling between the cavity and waveguides. 
We will focus on the situation where the frequency of 
the cavity lies between the band centers of the two 
waveguides, i.e. ujc ~ (oji + cj2)/2. In the numeri- 
cal calculation, we take wi = 9.5GHz, uji — 10.5GHz, 
6 = 6 = C = 0.3GHz, Cic = = Cc, and = lOGHz. 
We should first compare the above analytical BM solu- 
tion in the weak coupling limit with the exact numerical 
solution of the cavity field amplitude, the cavity photon 
number and the photoncurrent, given by Eqs. (1631) and 
and Eqs. ^ and (gS]), respectively. The resuh is 
plotted in Fig. [3] where the coupling rate rj — S^/S^q — 0.5 
which belongs to a weak coupling and the BM approxima- 
tion is applicable as we have shown in 2J| . From Fig. [3l 
we see that in the weak coupling limit, the BM solution is 
in good agreement with the exact solution, in particular 
when the driving frequency a;^ = uji. When the driving 
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FIG. 3: Comparison of the analytical solution in the weak cou- 
pling limit with the exact numerical solution at the coupling 
rate rj = ^/^o ~ 0.5 with different driving frequency, (al)-(al) 
and (a2)-(b2) are the cavity field amplitude and cavity pho- 
ton number, respectively, with the BM solution (red dashed 
line) and the exact solution (blue solid line). (a3)-(b3) are 
the photocurrent Ii{t) and l2{t) with the exact solution (red 
solid line and blue dashed line) and the BM sohition (blue 
dotted line and red dotted-dashed line). 



frequency equals to the cavity mode frequency, the de- 
viation between the exact solution and the BM solution 
becomes relatively large in the steady limit. This can be 
seen from Eq. (j63p where the amplified cavity field ampli- 
tude in the BM limit. En = , ^" , becomes sensitive 

to the approximation. But the qualitative behavior of 
the BM solution is still in good agreement with the exact 
solution. 

Next, we shall present the exact numerical solutions 
for different values of the coupling constants between 
the cavity and the waveguides, to examine the differ- 
ent photonic transport dynamics in the weak coupling 
as well as in the strong coupling regime. Fig. 2] shows 
the cavity field amplitude of Eq. ([5^ with different cou- 
pling strengths and different driving field frequencies. 
Fig. and Fig. |4l[b) correspond to the cases of the 
driving frequency being in resonance with the band cen- 
ter of CROW 1 and the cavity mode, respectively. Dif- 
ferent curves in each plot correspond to different cou- 
pling strengths between the cavity and waveguides. In 
the weak coupling regime, the behavior of the cavity field 
amplitude highly depends on the driving field frequency. 
When the driving field is not in resonance to the cavity 
mode, the field amplitude oscillates (as a result of the 
superposition of the driving field with the cavity mode) 
and decay (due to the dissipation induced by the cou- 
pling of the cavity with waveguides). When the driving 
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FIG. 4: The time evolution of the cavity field amplitude with 
different coupling strengthes: 77 = 0.5 (solid green line), 1.0 
(dashed blue line) and 2.0 (dotted red line), (a) The driving 
field frequency equals to the band center of the CROW 1. 
(b) The driving field frequency is in resonance to the cavity 
mode. 



field is in resonance to the cavity mode (w^ = Wc), the 
field amplitude increases gradually without oscillation. 
These numerical results agree with the BM solution of 
Eq. (|63ap . as shown in Fig. |3l In the strong coupling 
regime, the field amplitude keep oscillating without de- 
cay in both cases. The absence of the damping (dissipa- 
tion) is totally due to the non-Markovian memory effect, 
as we have pointed out recently [2^ [23] • The compli- 
cated oscillating behavior (dotted red curve) in Fig. IDJa) 
is an interference effect between the driving field and the 
cavity mode with Wd i^c- Only in the resonance case 
{iOd = uJc)-, the cavity field becomes coherent with the 
external driving field, as shown by the dotted red curve 
in Fig.mjb). 

Fig. [5] shows the thermal fiuctuation v{t,t) and the 
average photon number n{t) in the driven cavity for 
different coupling strengthes, different driving field fre- 
quencies and different initial temperatures. According 
to Eq. (j59bl) . when the cavity is empty initially, the av- 
erage cavity photon number is fully determined by the 
driving-induced field |j/(t)P and the thermal fluctuation 
induced correlation v{t,t). As shown in Fig. [SJ v{t,t) 
almost equals to zero when T — 5mK. In this case, 
n{t) = |(a(i))p ~ \y{t)\^, namely, the cavity field is 
a pure coherent field induced by the external driving 
field. When the temperature increases, the contribution 
of v{t,t) is increased as well. At T = 5K, the contribu- 
tion from the thermal fluctuation is still very small in the 
weak coupling, in comparison with the contribution from 
the driving-induced field y(t). However, v{t,t) becomes 
comparable with y{t) in the strong coupling, as shows in 
Fig. [5l[a3)-(b3). With further increase of the tempera- 
ture, thermal fiuctuation would destroy the coherence of 
the cavity field if the driving field is weak. 

Fig. |6] shows the photocurrents Ii (t) and I2 (t) fiowing 
from the cavity into the CROWs in different driving fre- 
quencies and coupling strengthes. The frequency of the 
external driving field is the crucial factor for the control 
of the photocurrents. When the external driving field is 
in resonance to the cavity field, the photocurrents flow- 
ing through CROW 1 and CROW 2 are equal because 
of the symmetric configuration, as shown in Fig. |ni[a2)- 
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FIG. 5: The time evolution of the thermal fluctuation v{t,t) 
(red dashed and dotted-dashed curves) and the cavity inten- 
sity (blue solid and dotted curves) at different initial tem- 
perature (T = 5K and T = 5mK ) with different coupling 
strength. The left plots correspond to ujd — uji and the right 

plots for Ud ~ i^a- 



(c2). However, when the driving field is in resonance 
to the band center of CROW 1, the photocurrent I\(t) 
flowing through CROW 1 is dominated while the pho- 
tocurrent l^ify flowing through CROW 2 is suppressed, 
see Fig.[SJal)-(cl). Similarly, when the driving field is in 
resonance to the band center of CROW 2, the photocur- 
rent Ii (t) becomes dominated and the photocurrent Ii (t) 
is suppressed, as shown by Fig. [6lja3)-(c3). These results 
explicitly demonstrate the controllability of the photonic 
quantum transport through the driven cavity. 

The above numerical analysis show that the photonic 
dynamics and transport in an all-optical circuit can be 
manipulated efficiently through the external driving sig- 
nals and the internal cavity-waveguide couplings. Pho- 
ton dissipation and thermal fluctuation appear in all the 
cases. In the weak coupling limit, the cavity mode will 
decay so that only the driving-induced field plays the 
role in the photonic transport. The strong coupling be- 
tween the cavity and the waveguides can largely suppress 
the cavity photon dissipation, due to the non-Markovian 
memory effect. In this case, the cavity mode can serve 
as a new control field in the photonic transport. The 
photonic thermal fluctuation can be suppressed by ei- 
ther lowing the temperature or making the device work 
in the high frequency region. These general properties 
provide a theoretical basis for the further development of 
integrated quantum photonic circuits. 
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FIG. 6: Photocurrent Ii{t) (solid red line) and hit) (dashed blue line) in different coupling strength and driving frequencies. 
Here, the initial temperature T — 5K. 



VI. SUMMARY AND PROSPECTIVE 

In summary, we have established a quantum trans- 
port theory to describe photonic transport in photonic 
networks. The photonic networks consist of all-optical 
circuits incorporating photonic bandgap waveguides and 
driven resonators. The transient transport of photons 
in photonic networks is controlled through the driven 
resonators. The photonic dynamics of the driven res- 
onators is determined by the master equation which is 
derived by treating the waveguides as an environment. 
The back-reactions between waveguides and resonators 
and thereby the dissipation and fluctuation arisen from 
the back-reactions are fully taken into account. The pho- 
tocurrents flowing from the resonators into the waveg- 
uides that describe the photonic transport in the network 
are obtained directly from the master equation. The the- 
ory can be applied to photonic transport in photonic net- 
works involving many photons as well as single photon. 

Comparing with the electron transport in mesoscopic 
systems, we have shown several crucial differences for the 
photonic transport. Instead of the bias and gate voltage 
controls in electron transport, the photonic transport is 
manipulated through external driving fields applying to 
waveguides and resonators. Differing from the electron 
transport where the bias and gate voltages controls are 



imbedded into the spectral densities or the energy lev- 
els of the system, the external driving fields applying to 
waveguides and resonators turn out to be the modified 
driving fields acting explicitly on the resonators, as shown 
in the master equation of Eq. (|27p . The resulting lesser 
Green function, see Eq. (^5]) . and thereby the photocur- 
rent, see Eq. (|49p . contain additional terms induced from 
the external driving fields that explicitly determine the 
photonic transport in the photonic network. In other 
words, the present theory shows that the transport con- 
trols through the external fields behave very differently 
for electrons and photons. Furthermore, the explicit ini- 
tial state dependence also makes the present theory be- 
come more powerful in the study of transient transport 
phenomena. 

As a simple illustration of the theory, we investigate 
the photonic transport dynamics of the driven nanocav- 
ity coupled to two waveguides. The cavity field, the 
cavity photon number as well as the photocurrents flow- 
ing through the waveguides are solved analytically in the 
weak coupling limit and are also exactly calculated nu- 
merically for strong couplings with different driving fre- 
quencies. Non-Markovian memory effects in the strong 
coupling regime are shown explicitly. The thermal fluctu- 
ation effect to the coherent property of the cavity field is 
also demonstrated. Moreover, the analytical and numer- 
ical results show the strong controllability of wavelength 
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selective transport. The controllability together with the 
simplicity of the photonic circuit implementation in pho- 
tonic crystals make it useful in the further development 
of integrated quantum photonic circuits. 

We should further point out that the present the- 
ory can also be directly extended to investigate vari- 
ous photonic transport phenomena in other electromag- 
netic matematerials 160,1 and phonon transport in heat- 
conducting systems [6l[. The generality of the present 
theory incorporating with the Feynman- Vernon influence 
functional approach and the Keldysh's nonequilibirum 
Green function technique together also makes it pow- 
erful in the investigation of open quantum system, in 
particular the nonequilibrium dynamics. The generaliza- 
tion of the present theory to the nonequilibirum dynam- 
ics of ultracold atomic Boson-Einstein condensation is in 
progress. Various novel devices, such as photon entangle- 
ment through photonic crystal waveguides |6^ l63j , can 
be designed with the help of the present theory. And 
more applications will be presented in future works. 
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Appendix A: Derivation of the propagating function 



The influence functional of Eq. ([2T|) modifies the 
original action of the system into an effective one, 

gi(Ss[a*,a]-S's[Q:'*,a'])^jQ,*Q,. Q,/*Q,/J _ ^j-S,{t[a' a;a" a'] 

which dramatically changes the dynamics of the driven 
resonators. The explicit change is manifested through 
the generating function of Eq. (I19[) by carrying out 
the path integral with respect to the effective action 
Ses[ot.* ot', (y.'*a']. While the path integral V[a.*a.; a.'*cx'] 
integrates over all the forward paths q:*(t), q:(t) and the 
backward paths Q.'*{T),a'{T) in the Bergman complex 
space bounded by a*{t) — a*p a.{to) = ccq and a.'* (to) = 
a'l,a.'{t) = a'p respectively. Since Scs[ot*a.;a'*a'] is 
only a quadratic function in terms of the integral vari- 
ables, the path integrals of Eq. (|19l) can be reduced to 
Gaussian integrals so that we can use the stationary 
path method to exactly carry them out [gj] . The result- 
ing propagating function is a function of the stationary 
paths: 

J{aj,a'^,t\at,a'*,to) = ^(i) exp{i[aJ.Q;(i) +oi''{tQ)a, 



+ /t(r)(a'(r)-a(r))]} , 



[K^(r)-«t(r))/(r) 



where A{t) is the contribution arisen from the fluctu- 
ations around the stationary paths and is given after 
Eq. (f24| . The stationary paths obey the equations of 
motion : 



da{T) 
dT 



+ iu}{T)a{T)+ f dT'g{T,T')a{T') 

J to 

+ f dT'g{T,T')[a{T') - c'(t')] ^ -tfir), (A2a) 

Jto 



da'ir) 
dT 



- iu}{T)a.' (t) 



- /*dr'g(r,T')a'(T')+ f dT'^{Ty )cx{T') 

Jt Jto 

+ rdT'g(r,r')[a(T')-a'(r')] = -z/(r), (A2b) 

Jto 

subjected with the boundary conditions Q:(to) = OLi and 
a.'{t) — a'y. The equations of motion for ot^{T) and 

ol'\t) follow by exchanging a. and a' in Eq. (jA2p and 
taking a conjugate transpose, subjected to the boundary 
conditions a'^(io) = ol'\ and ot^it) = q:|. 

To express the master equation independent of the co- 
herent state representation, we shall further factorize the 
boundary values of the stationary paths, a(to) = ccq, 
a.'{t) — a.'j, through the following transformation: 



a'(T) — q:(t) = u{T,t){a.'f — Q.{t)), 



(A3a) 



to 



(Al) 



a(T) = m(t, to)aQ + v{t, t){a'j: - a{t)) + y{T) (A3b) 

and a similar transformation for conjugate variables 
(with the exchange of a. with ol' for the boundary values 
a.^{t) = ol'j and a'^(to) = ol'\), where M(r, io), u{T,t), 
v{T,t) are by A^ matrices, and ^(t) is 1 x A^ matrix, 
A^ is the dimension of the cavity system. Substituting 
the solutions of Eq. (jA3[) into the equations of motion of 
Eq. (|A2[) and comparing with the boundary values, we 
obtain the equations of motion for u(t, to), u(t, t), v(r, t) 
and y(r), given by Eq. (^5]) . By multiplying u'^{T,to) to 
Eq. (|25ap and integrating r from to to t, we obtain the re- 
lation u{t,to) = u^{to,t). Similarly, we can also directly 
obtain the solutions of Eqs. (|25cp and (j25dp . which are 
given by Eq. (^5]) . 

Now, let T = to in Eq. (|A3ap and t = i in Eq. (jA3bp . 
we can express a{t) and a' (to) in terms of the boundary 
conditions cto and a'^: 

cx{t) =w{t)[u{t,tQ)ai + v{t,t)a'f +y{t)] , (A4a) 
a' (to) =u^{t, to)[l - w{t)v{t, t)]a^. - u^t, to)w{t)y{t) 
+ [l-u^{t,to)w{t)u{t,to)]ao , (A4b) 

here w{t) = [1 + v{t,t)]-'^ = w^t). Similarly, a^to) 
and ot'\t) can be obtained by exchanging at and a' in 
Eq. (|A4[) and taking a conjugate transpose. Substituting 
these results with Eq. (jA3p into Eq. (jAip . we obtain the 
final form of the propagating function for the reduced 
density matrix given by Eq. 
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Appendix B: Analytical solution in weak coupling 
limit with external driving field 

In this appendix, we shall present an analytical solu- 
tions in the weak coupling limit for the photonic net- 
work concerned in Sec. V. First, let us solve the photon 
propagating function u{t,tQ). To do so, let u(t,ta) = 
e-^"-(*-*o)A(i). Then Eq. d^^a]) is reduced to 

r dt'g(t')e'^^''A(t~t') = 0. (Bl) 
at Jo 

When the spectrum of the waveguides are broad enough 
and the coupling between the waveguides and the sys- 
tem is weak, the memory effect between the system and 
the reservoir can be ignored so that the Markov limit is 
reached. In other words, A(t — r) in the integration can 
be replaced by A(t) and the integration time t' can be 
considered much shorter in comparison with the charac- 
ter time t of the cavity field 47] . Then, 



f-to 



dt'^{t')e"^''' A{t-t') 



lim 



t-to 



(B2) 

X]a=i is the cavity frequency shift 



r 57^^- The cavity damping rate 



where Sujc 
with 5uja - 

K = X]q=i with Ka — Ja{uJc)/2. Thus, the photonic 
propagating function is approximated analytically by 

M(i,to) e-(''"^+'')(*-*«\ (B3) 

and the renormalized cavity frequency uj'^ = ujc + 5ujc- 
Next, we shall solve the correlation function v{t,t) in 

the same approximation. Note that 7^ ^^^^|r=t, 

one should not directly use Eq. (|25cp to solve v{t, t). It is 
more convenient to take the time derivative of the solu- 
tion v{t, t) given by (|26cp and utilize the Eq. (|25ap . which 
leads to 

2Re / dTg{t~ T)v{T,t) 

J to 

= 2Re / dT^{t ~ t)u{t, t). (B4) 

J to 

Then using the same approximation of Eq. (IB2p . we have 



dt 



Re / dTg{t-T)v{T,t) 



to 



Re / dTg{t — T)u{T,t) 



Re / dTg{t 

J to 
Kv{t, t), 

Re / dTg{t 

J to 



re 



v{t,t) 



c^Y^JMnM. (B5) 

a=l 

Thus, the solution of Eq. (jB4p is approximately given by 
«(t,i)^nK,r)[l-e-2'=(*-*«)], (B6) 



where n{oj'^,T) = "'ct('^c)'^a('^c)/2'*- The analyti- 

cal solutions of Eqs. (|B3p and (|B6p are the well-known 
Born-Markov limit of the cavity field coupled to a ther- 
mal bath. 

With the external driving field f{t) = ^Qg"*'^'***-*"' 
being explicitly added, the equation of motion for the 
driving-induced field j/(t) becomes: 



dyjt) 

dt 



+ iujMt) + / dTg(t - T)y{T) = -iEoe~''^'''^*-*°\ 

J to 

(B7) 

The homogeneous solution of this equation is just u(r, to) 
with the character frequency uJc, whose BM limit solution 
has been given by Eq. (jB3p . The inhomogeneous solution 
of y{t) must have a form _Be~*'^'*(*~*°\ where S is a con- 
stant and ujd is the character frequency. In the BM limit, 
using the same approximation but noting the different 
character frequency, we find that 



where 



B 



tan 



Eo cxp (— i(/ 



(B8) 



T,a = l ^'^o 



with the driving field induced frequency shift SCja — 
^/o°" and k = ELi with = UuJ,)/2. 

Put the homogeneous solution and the inhomogeneous 
solution together with the initial condition y{to) = 0, we 
obtain the driving-induced field in the BM limit: 

(B9) 

From the above solution for the propagating function 
m(t, io), the correlation function v{t,t) and the driving- 
induced field y{t), it is easy to find analytically the cav- 
ity field and the cavity photon number in the BM limit, 
i.e. Eq. (|63p . Utilizing the similar approximation in the 
above derivation, the photocurrent flowing through the 
waveguide a = 1, 2 in the BM limit are given by Eq. (|64p . 
Furthermore, we can also find the source term in Eq. p4)) 
in the BM limit by substituting the BM solution of the 
cavity field, Eq. (|55al) . into Eq. (jggaP : 



S(t) 



2kEl 



! ~ .2^— {l + e--(*-*«)cos[(^,-c.^)(i-io)] 
+ ii^^i^e-«(*-*«) sin [{ud - u',){t - to)]]. 

(BIO) 



As a self-consistent check, the BM solutions of the cavity 
photon number, the driving source and the photocurrent 
fiowing into waveguides satisfy the quantum continuous 
equation: 



dn{t) 
dt 



s{t)-j2ut). 



(Bll) 
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